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Summary 

Bats and rodents are being increasingly recognized as reservoirs of emerging zoonotic 
viruses. Various studies have investigated bat viruses in tropical regions, but to date 
there are no data regarding viruses with zoonotic potential that circulate in bat and rat 
populations in Viet Nam. To address this paucity of data, we sampled three bat farms 
and three wet markets trading in rat meat in the Mekong Delta region of southern Viet 
Nam. Faecal and urine samples were screened for the presence of RNA from para- 
myxoviruses, coronaviruses and filoviruses. Paramyxovirus RNA was detected in 4 of 
248 (1%) and 11 of 222 (4.9%) bat faecal and urine samples, respectively. Coronavirus 
RNA was detected in 55 of 248 (22%) of bat faecal samples; filovirus RNA was not 
detected in any of the bat samples. Further, coronavirus RNA was detected in 12 of 
270 (4.4%) of rat faecal samples; all samples tested negative for paramyxovirus. 
Phylogenetic analysis revealed that the bat paramyxoviruses and bat and rat coronavi- 
ruses were related to viruses circulating in bat and rodent populations globally, but 
showed no cross-species mixing of viruses between bat and rat populations within 
Viet Nam. Our study shows that potentially novel variants of paramyxoviruses and 
coronaviruses commonly circulate in bat and rat populations in Viet Nam. Further 
characterization of the viruses and additional human and animal surveillance is re- 
quired to evaluate the likelihood of viral spillover and to assess whether these viruses 


pose a risk to human health. 
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1 | INTRODUCTION 


A large proportion of the agents of emergent infectious disease have 
a zoonotic origin. These zoonotic pathogens fall into a wide spec- 
trum of genera, but due to their high genetic variability and wide cir- 
culation, RNA viruses arguably pose the most significant threat to 
human health. Further, due to the close proximity between animals 
and humans in low- and middle-income countries, human popula- 
tions in these locations are disproportionally at risk of exposure from 
these viral pathogens (Paterson et al., 2014). Although animal viruses 
with zoonotic potential have likely been circulating continually, a 
high number of zoonotic RNA viruses have been discovered in recent 
years. This phenomenon is not only dependent on more thorough 
and enhanced detection methods but is also likely associated with 
the changing behaviour of human populations and the closer proxim- 
ity between humans and the animals that act as reservoirs for these 
viruses (Calisher, Childs, Field, Holmes, & Schountz, 2006; Moratelli 
& Calisher, 2015; Sasaki et al., 2012; Wong, Lau, Woo, & Yuen, 2007). 

Multiple studies have implicated bats to be the most likely reser- 
voir of numerous zoonotic viruses (Calisher et al., 2006; Moratelli & 
Calisher, 2015; Wong et al., 2007), and there are a range of contrib- 
uting factors that may make so-called “spillover” events, in which hu- 
mans or other animals are infected, more probable. Many bat species, 
including the insectivorous Scotophilus kuhlii present across South and 
South-East Asia (Bates et al., 2016), live in large groups with a range 
of social behaviours involving close and prolonged contact with others 
in the roost; this activity may facilitate the horizontal transmission of 
viruses between roost members. In addition, some bat species have a 
long life span (up to 35 years) and are able to travel long distances (e.g., 
Eidolon helvum can migrate >4,500 km), thus increasing the likelihood 
of exposure to infectious agents (Drexler, Corman, & Drosten, 2014; 
Escaffre, Borisevich, & Rockx, 2013; Plowright et al., 2015). The mod- 
ification and destruction of natural habitats increase the likelihood of 
contact between bats and humans, thus providing new opportunities 
for interspecies viral exchange. 

Bats are a known reservoir for rabies virus (Escobar et al., 2015; 
Jakava-Viljanen et al., 2015; Moratelli & Calisher, 2015; Rocha, de 
Oliveira, Heinemann, & Goncalves, 2015), and have recently been 
associated with other viral pathogens with a severe infection phe- 
notype in humans. Notoriously, Nipah and Menangle viruses, both 
paramyxoviruses, have caused outbreaks in humans, horses and pigs 
in Australia, Cambodia, Malaysia, Bangladesh and India (Escaffre et al., 
2013; Kulkarni, Tosh, Venkatesh, & Senthil, 2013). Furthermore, SARS 
coronavirus (SARS-CoV) has been identified in Chinese bats, and it has 
been shown that the SARS-CoV genome sequence generated from 
humans and civets during the 2002-2003 outbreak in China phyloge- 
netically clustered within the bat associated group of SARS-CoV-like 
viruses (Calisher et al., 2006; Drexler et al., 2014; Moratelli & Calisher, 
2015; Wong etal., 2007). Filoviruses, including Ebola and Marburg 
viruses, also pose a significant threat to human health, and despite 
human outbreaks being rare and sporadic, they are associated with 
a high case fatality rate (Hoffmann et al., 2013; Jayme et al., 2015; 
Kgaladi et al., 2013; Plowright et al., 2015). 


Impacts 

e Bats and rodents are known reservoirs of highly diverse 
viral and bacterial populations, and a number of these vi- 
ruses have been implicated in the emergence of novel 
infectious diseases in humans. 

e The close proximity of humans to bat and rodent popula- 
tions in South-East Asia creates frequent opportunities 
for viral spillover and thus poses an unknown risk to 
human health. 

e Viral surveillance in animal reservoirs is an important step 
to understanding the exposure of humans to potential 
zoonoses, the types of human-animal interaction that im- 
pact potential for spillover infection and the factors that 
determine the transmissibility and pathogenicity of viral 


zoonoses in humans. 


Viet Nam is a middle-income country in South-East Asia. Bats are 
common in Viet Nam and may pose a threat from the circulation of 
zoonotic RNA viruses. Further, in parts of the country, rats are com- 
monly trapped in the rice fields and sold live for consumption; 3,300- 
3,600 tonnes are sold nationally annually (Van Cuong et al., 2015). 
Rats and other rodents are also recognized as a reservoir of zoonotic 
viruses that can be transmitted via close contact with saliva, urine or 
faeces; the circulation of hantavirus has been reported in rats in Viet 
Nam (Van Cuong et al., 2015). As South-East Asia is a hot spot for zoo- 
notic viruses, the WT-VIZIONS (Wellcome Trust-Vietnamese Initiative 
on Zoonotic Infections) project (Rabaa et al., 2015) is aiming to gen- 
erate data on the circulation of viral zoonotic pathogens that pose a 
risk to human health. With bats suspected to be the main reservoir 
for more than 200 viral species, 248 bat samples from three differ- 
ent guano farms in the south of Viet Nam (Dong Thap Province) were 
screened for paramyxoviruses, coronaviruses and filoviruses to assess 
the potential risk to human health. In addition, 270 rodent faecal sam- 
ples from the same province were screened for paramyxoviruses and 


coronaviruses. 


2 | MATERIALS AND METHODS 


2.1 | Bat samples 


Insectivore bat colonies (S. kuhlii) rather than individual bats were tar- 
geted for screening; three bat guano farms were identified in three 
geographically distinct sampling locations in the south of Viet Nam. 
These farms were selected because they were the only bat guano 
farms identified and consenting to participate in regular sampling 
within the catchment area for Dong Thap Provincial Hospital, a pri- 
mary site of ongoing human sampling for the identification of zoonotic 
infections under the VIZIONS project (Rabaa et al., 2015). The three 
sampling locations were visited every 12 weeks (one sampling event) 


on six occasions between 2012 and 2014. The farms were visited in 
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the afternoon; to catch urine and faeces, a ten-square metre section 
of plastic sheeting was placed under trees in which bats were roosting. 
The urine and faecal samples were collected from the plastic sheet- 
ing using a sterile Pasteur pipette for the urine and a plastic scoop 
for the faeces. These samples were then placed in a graduated 2-ml 
tube (Sarstedt microtubes, Denmark) and transferred to the labora- 
tory in Ho Chi Minh City for processing and/or storage at -80°C. At 
each sampling event, we aimed to collect 10-15 samples per farm. 
Depending on the roosting density and activity of the bats, we esti- 
mated the plastic sheeting could hold urine and faeces for up to 200 


individual bats. 


2.2 | Rat samples 


The rodent surveys and sampling were conducted as previously de- 
scribed by Van Cuong et al. (2015). Briefly, rats were purchased from 
three markets in Dong Thap Province in the Mekong Delta (two rat 
markets in Cao Lanh City and one market in Tam Nong District) every 
4 months for 2 years, between 2012 and 2014. Specimens used in this 
investigation included faeces from 270 individual rats of four differ- 
ent rat species: Bandicota indica (n = 8), Rattus argentiventer (n = 234), 


Rattus losea (n = 20) and Rattus tanezumi (n = 8). 


2.3. | RNAextraction, PCR screening and sequencing 


Total nucleic acid was extracted from the bat urine, bat faecal and rat 
faecal samples using an automated extracted system (MagNA Pure 
96 System, Roche) with MagNA Pure 96 Viral RNA Small Volume Kits 
(Roche). Following extraction, nucleic acid samples were screened for 
target viruses using pan-family primers for paramyxoviruses, coro- 
naviruses and filoviruses. The primer sequences of viral targets and 
the nested PCR conditions were as previously described (Heaton 
et al., 1997; Ogawa et al., 2011; Poon & Peiris, 2008; Tong, Chern, 
Li, Pallansch, & Anderson, 2008). Positive controls were included for 
each run and provided by CSIRO Biosecurity Flagship & Australian 
Animal Health Laboratory, Geelong, Victoria, Australia. All PCR 
amplification-positive samples were subsequently characterized by 
Sanger sequencing using an ABI 3700 sequencer in both forward and 
reverse directions using the primers from the second round nested 


PCR of each target virus. 


2.4 | Phylogenetic analysis 


For both paramyxovirus and coronavirus, partial RdRp sequences were 
aligned with subsets of publicly available reference sequences, yielding 
an 84-sequence data set for paramyxoviruses (344 bp; sites 14210- 
14553 within the genome), a 153-sequence data set for alphacoro- 
naviruses and a 93-sequence data set for betacoronaviruses (407 bp; 
sites 15149-15565 within the genome), using Seqotron (Fourment 
& Holmes, 2016). Identical sequences and those with greater than 
1% ambiguity or missing sequence data were subsequently removed 
from the alignments. jModelTest was run, and indicated the GTR+I+G 


model to be the best-fit model of nucleotide substitution for all three 
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data sets (Posada, 2009). Maximum-likelihood phylogenies were run 
using RAxML under the GTRGAMMAI model with 1,000 bootstrap 
replications for each dataset (Stamatakis, 2006). Trees were visualized 


and annotated using FigTree (v1.4.2). 


3 | RESULTS 


3.1 | Screening of bat samples for paramyxovirus, 
coronavirus and filovirus 


Bat urine and faeces collected from three guano farms in the south 
of Viet Nam were extracted and subjected to PCR amplification for 
paramyxoviruses, coronaviruses and filoviruses. Four of 248 (1.6%) bat 
faecal samples and 11 of 222 (4.9%) urine samples tested positive for 
paramyxovirus RNA. The paramyxovirus-positive faecal samples were 
identified in sampling round two ona single farm and in sampling round 
three on two farms. The positive urine samples were detected on three 
bat guano farms during the sixth round of sampling. Screening of bat 
faeces and urine for coronaviruses indicated that 55 of 248 (22.2%) of 
faecal samples were positive for coronavirus RNA; none of the urine 
samples contained detectable coronavirus RNA. Coronavirus RNA was 
detected in bat faecal samples isolated from all three sampled farms 
at all of the sampling time points. The bat urine and faecal samples 
were also screened for filoviruses, but no PCR amplification for filo- 
virus RNA was detected despite amplification of appropriate controls. 
Specific information on coronavirus and paramyxovirus positivity by 


sampling site, time and species is indicated in Tables S1-S3. 


3.2 | Screening of rat samples for coronavirus and 
paramyxovirus 


Of the 270 rat faecal samples screened for coronaviruses and para- 
myxoviruses, 12 of 270 (4.4%) were positive for coronavirus RNA, 
while none of the samples had detectable paramyxovirus RNA. The 
cDNA amplicons were sequenced and all twelve positive samples 
produced sequences suitable for phylogenetic analysis. Notably, all 
coronavirus-positive faecal samples were obtained from rats of one 


species, R. argentiventer. 


3.3 | Phylogenetic relationships among Vietnamese 
coronaviruses and paramyxoviruses from bats and rats 


Fifteen RNA sequences obtained from the bat paramyxovirus RNA- 
positive samples (4/248 faecal samples and 11/222 urine samples) 
(GenBank accession numbers from KX092148 to KX092159; Table S4) 
were aligned with 69 reference sequences from the conserved region 
within the RNA-dependent RNA polymerase (RdRp) gene (344 bp) 
accessed from GenBank. A maximum-likelihood (ML) phylogeny sug- 
gested that all of these viruses were closely related and clustered with 
paramyxoviruses found in bat populations sampled globally, although 
low bootstrap values across much of the phylogeny indicate a largely 
unresolved evolutionary history of these viruses (Figure 1) using a 


short fragment of the RdRp gene alone. Based on ICTV classification, 
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FIGURE 1. The phylogenetic relationships of paramyxoviruses sampled from Vietnamese bats. Maximum-likelihood phylogeny constructed 
using RNA sequences from 15 Vietnamese bat paramyxovirus RNA-positive samples and 69 reference sequences from the conserved region 


within the RNA-dependent RNA polymerase (RdRp) gene (344 bp). The tree is mid-point-rooted for the purpose of clarity. Scale bar indicates the 
number of substitutions per site. Bootstrap support values are shown for nodes with 280% bootstrap support. Sequences from Vietnamese bats 


are indicated in red 
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TABLE 1 Maximum nucleotide and amino acid identities of Vietnamese bat paramyxovirus sequences relative to other known 


paramyxoviruses 


Highest nucleotide 


Highest amino 


Paramyxovirus identity (%) Closest match acid identity (%) Closest match 
05VZ_75_55_L13_R6_VN 80.18 KC578579|Kerivoula argentata|South 92.11 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_55_L09_R3_VN 78.91 KC578575|Nycteris thebaica|South 93.75 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_55_L10_R6_VN 80.47 KC578572|Eptesicus 92.94 KC578572|Eptesicus 
hottentotus|South Africa hottentotus|South Africa 
05VZ_75_98_LO2_R3_VN 79.26 KC578575|Nycteris thebaica|South 93.86 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_LO05_R6_VN 80.78 KC578579|Kerivoula argentata|South 93.14 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_L06_R2_VN 79.26 KC578575|Nycteris thebaica|South 93.86 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_LO8_R6_VN 80.82 KC578579|Kerivoula argentata|South 92.45 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_LO9_R6_VN T33 KC578572|Eptesicus 91.23 KC578572|Eptesicus 
hottentotus|South Africa hottentotus|South Africa 
05VZ_75_99_L13_R6_VN 80.97 KC578579|Kerivoula argentata|South 91.82 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_L10_R6_VN 77.78 KC578575|Nycteris thebaica|South 90.72 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_98_L15_R6_VN 80.93 KC578572|Eptesicus 94.64 KC578572|Eptesicus 
hottentotus|South Africa hottentotus|South Africa 
05VZ_75_99_L01_R6_VN 80.68 KC578572|Eptesicus 92.75 KC578572|Eptesicus 
hottentotus|South Africa hottentotus|South Africa 
05VZ_75_99_LO2_R6_VN 80.84 KC578575|Nycteris thebaica|South 94.06 KC578572|Eptesicus 
Africa hottentotus|South Africa 
05VZ_75_99_LO05_R6_VN 81.97 KC578579|Kerivoula argentata|South 92.86 KC578572|Eptesicus 
Africa hottentotus|South Africa 


the bat paramyxoviruses determined in this study represented border- 
line novel species, with maximum nucleotide identities ranging from 
77.8% to 81.97% compared to publicly available sequences (Table 1). 
All of our sequences showed the highest identity to recently described 
viruses isolated from various bat species in sub-Saharan Africa. 
Sequence data indicated that all bat coronaviruses characterized 
in this study were of the alphacoronavirus lineage, while all rat coro- 
naviruses were of the betacoronavirus lineage. The alphacoronavirus 
sequences obtained from bat samples (n = 55; Table 2 and S5) and 
betacoronavirus sequences from rat samples (n = 12; Table 3 and S6) 
(GenBank accession numbers: KX092163-KX092228) were aligned 
with reference sequences (407 bp), identical sequences and lower 
quality sequences were removed, and ML trees were constructed. 
Based on the 2009 ITCV classification, none of the coronaviruses 
detected in this study represented novel species (Tables 2 and 3). 
Phylogenetic reconstruction indicated that bat alphacoronaviruses 
from Vietnam clustered with other bat viruses detected in China and 
the Philippines (Figure 2a), and showed the highest amino acid identity 
to bat alphacoronaviruses described in China (Table 2). The betacoro- 
naviruses detected in rats in this study clustered with reference sam- 


ples originating from other mammal populations, with one of the two 


clusters found to be closely related to betacoronaviruses detected 
in rodent populations sampled in China (Figure 3). These sequences 
showed the highest identity to rat and mouse betacoronaviruses from 
China (Table 3). Again, inference is somewhat limited by the length of 
these sequences and additional sequencing effort is needed to con- 


firm these preliminary findings. 


4 | DISCUSSION 


Here, we sought to investigate the prevalence of paramyxovirus, 
coronavirus and filovirus in bat and rat populations in southern Viet 
Nam. For this study, we implemented a sampling and testing protocol 
that sought to maximize the likelihood of detection if the target vi- 
ruses were present. Firstly, we collected pooled urine and faecal sam- 
ples from three different bat guano farms in the south of Viet Nam, 
rather than catching and sampling individual animals, so each sample 
potentially represented multiple animals. Bats forage nocturnally and 
roost communally in trees during the day; during this period, they fre- 
quently urinate and defecate, allowing collection on plastic sheeting. 


This approach has been used successfully (Field et al., 2015) and has 
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TABLE 2 Maximum nucleotide and amino acid identities of Vietnamese bat alphacoronavirus sequences relative to other known 


alphacoronaviruses 


Alphacoronavirus 


75_55_L01_R1_CoV_BAT_VN 


75_55_LO3_R3_CoV_BAT_VN 


75_55_LO06_R2_CoV_BAT_VN 


75_55_LO9_R6_CoV_BAT_VN 


75_55_L10_R3_CoV_BAT_VN 


75_55_L10_R6_CoV_BAT_VN 


75_55_L12_R6_CoV_BAT_VN 


75_55_L13_R6_CoV_BAT_VN 


75_55_L15_R2_CoV_BAT_VN 


75_98_LO7_R5_CoV_BAT_VN 


75_98_L01_R3_CoV_BAT_VN 


75_98_L05_R2_CoV_BAT_VN 


75_98_LO7_R2_CoV_BAT_VN 


75_99_L06_R4_CoV_BAT_VN 


75_55_L01_R4_CoV_BAT_VN 


75_55_LO2_R1_CoV_BAT_VN 


75_55_LO2_R2_CoV_BAT_VN 


75_55_LO3_R5_CoV_BAT_VN 


75_55_LO3_R6_CoV_BAT_VN 


75_55_L04_R1_CoV_BAT_VN 


75_55_L04_R5_CoV_BAT_VN 


75_55_LO5_R5_CoV_BAT_VN 


75_55_LO5_R6_CoV_BAT_VN 


75_55_LO9_R3_CoV_BAT_VN 


Highest 


nucleotide 


identity (%) 


95.55 


95.81 


95.81 


96.34 


96.07 


Ooi 


96.07 


96.34 


96.07 


94.85 


95.81 


Grol 


96.07 


ES) 


84.44 


96.30 


97.04 


PUS.Sib 


95.80 


96.30 


96.30 


PENS 


97.04 


84.69 


Closest match 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


AB539080|Scotophilus 


kuhlii 


Philippines 


DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 
DQ648821|Unknown bat 
species|China 


DQ648821|Unknown bat 
species|China 


Highest amino 


acid identity 
(%) Closest match 
97.78 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
97.78 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
97.44 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
97.78 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
98.52 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
93.33 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
99.26 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
100 NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
OP8}858) NC_009657|Scotophilus spp.|China, 
DQ648821|Unknown bat species|China 
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TABLE 2 (Continued) 


Highest Highest amino 
nucleotide acid identity 
Alphacoronavirus identity (%) Closest match (%) Closest match 
75_55_LO9_R5_CoV_BAT_VN 96.79 DQ648821|Unknown bat 99.26 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L11_R6_CoV_BAT_VN 96.78 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L13_R2_CoV_BAT_VN 85.68 DQ648821|Unknown bat 94.82 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L13_R4_CoV_BAT_VN 97.04 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L14_R4_CoV_BAT_VN 96.34 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L15_R3_CoV_BAT_VN 96.30 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_55_L15_R5_CoV_BAT_VN 96.79 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_LO5_R6_CoV_BAT_VN 84.44 DQ648821|Unknown bat 94.07 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_LO07_R4_CoV_BAT_VN 96.30 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L10_R4_CoV_BAT_VN 97.04 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L13_R4_CoV_BAT_VN 84.20 DQ648821|Unknown bat 91.85 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L14_R2_CoV_BAT_VN 97.04 DQ648821|Unknown bat 100.00 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L14_R4_CoV_BAT_VN 96.54 DQ648821|Unknown bat 97.78 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L15_R4_CoV_BAT_VN is) Chil DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_98_L15_R5_CoV_BAT_VN 96.40 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_LO2_R3_CoV_BAT_VN 95.56 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L03_R1_CoV_BAT_VN 96.79 DQ648821|Unknown bat 99.26 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_LO06_R1_CoV_BAT_VN 97.04 DQ648821|Unknown bat 9926 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_LO7_R5_CoV_BAT_VN 96.66 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_LO9_R1_CoV_BAT_VN 96.30 DQ648821|Unknown bat CN) XS) NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L12_R3_CoV_BAT_VN 95.56 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L14_R1_CoV_BAT_VN 96.30 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L15_R3_CoV_BAT_VN 95.56 DQ648821|Unknown bat 100 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L15_R5_CoV_BAT_VN 96.27 DQ648821|Unknown bat O95 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
75_99_L01_R5_CoV_BAT_VN 84.44 DQ648821|Unknown bat 94.07 NC_009657|Scotophilus spp.|China, 
species|China DQ648821|Unknown bat species|China 
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TABLE 3 Maximum nucleotide and amino acid identities of Vietnamese rat betacoronavirus sequences relative to other known 


betacoronaviruses 


Highest nucleotide 


Highest amino acid 


Betacoronavirus identity (%) Closest match identity (%) Closest match 
75_62_L01_R6_RAT_CoV_VN 97.54 KF294357|Apodemus 100 KF294357|Apodemus agrarius|China, 
agrarius|China NC_026011|Rattus 
norvegicus|China 
75_62_L12_R2_RAT_CoV_VN 93.86 KF294372|Niviventer 98.52 KF294370|Rattus tanezumi|China, 
confucianus|China KF294372|Niviventer 
confucianus|China 
75_63_L04_R5_RAT_CoV_VN 93.86 KF294372|Niviventer 98.52 KF294370|Rattus tanezumi|China, 
confucianus|China KF294372|Niviventer 
confucianus|China 
75_63_L11_R6_RAT_CoV_VN 90.66 KF294372|Niviventer 94.82 KF294370|Rattus tanezumi|China, 
confucianus|China KF294372|Niviventer 
confucianus|China 
75_65_LO2_R6_RAT_CoV_VN 97.54 KF294357|Apodemus 99.26 KF294357|Apodemus agrarius|China, 
agrarius|China NC_026011|Rattus 
norvegicus|China 
75_65_LO3_R6_RAT_CoV_VN 30) KF294372|Niviventer 97.04 KF294370|Rattus tanezumi|China, 
confucianus|China KF294372|Niviventer 
confucianus|China 
75_65_LO7_R5_RAT_CoV_VN 97.79 KF294357|Apodemus 99.04 KF294357|Apodemus agrarius|China, 
agrarius|China NC_026011|Rattus 
norvegicus|China 
75_65_LO7_R6_RAT_CoV_VN 97.00 KF294357|Apodemus 98.52 KF294357|Apodemus agrarius|China, 
agrarius|China NC_026011|Rattus 
norvegicus|China 
75_65_LO9_R2_RAT_CoV_VN 93.12 KF294372|Niviventer 96.30 KF294370|Rattus tanezumi|China, 
confucianus|China KF294372|Niviventer 
confucianus|China 
75_65_L10_R5_RAT_CoV_VN Serge KF294357|Apodemus 100 KF294357|Apodemus agrarius|China, 


agrarius|China 


been shown to increase the number of animals sampled and there- 
fore increase the chance of detecting viral nucleic acid. A limitation 
of this investigation is that the pooling of samples from many animals 
may have overestimated the true prevalence, given the possibility that 
multiple positive (pooled) samples may reflect a single individual. 

The methods used in the current study resulted in the identi- 
fication of previously unidentified paramyxoviruses and coronavi- 
ruses in Viet Nam and has expanded our knowledge of circulation of 
viruses with zoonotic potential in the bat and rat population in Viet 
Nam. Recently, bats have been shown to harbour a great diversity of 
previously unknown paramyxoviruses, some of which have been as- 
sociated with zoonotic events (Baker et al., 2013). The interrelation- 
ship of these paramyxoviruses with bats highlights the complex role 


that bats play as a reservoir for these viruses. Despite our sampling 


NC_026011|Rattus 
norvegicus|China 


framework being narrow, it was relatively thorough, as we screened 
248 bat faecal samples and 222 bat urine samples. RNA from novel 
paramyxoviruses was detected by PCR and could be characterized by 
sequencing in approximately 2% of bat faecal samples and nearly 5% 
of bat urine samples. Our study and studies by others indicate that it 
is highly likely that there are more extant bat-derived paramyxoviruses 
yet to be characterized (Drexler et al., 2009, 2012; Lau et al., 2010). 
For this reason, it is crucial that the potential consequences and risk of 
zoonotic spillover of such viruses are investigated further. 

Our data additionally suggested that alphacoronaviruses are 
constantly circulating in the bat population in Vietnam. The pres- 
ence of coronavirus RNA in 22% of bat samples showed that 
coronaviruses may be endemic, with highly similar viruses circulat- 


ing on all of the sampled farms at different time points. Ge et al., 


FIGURE 2 The phylogenetic relationships of alphacoronaviruses sampled from Vietnamese bats. Maximum-likelihood phylogeny constructed 
using RNA sequences from 40 Vietnamese bats and 98 reference sequences from the conserved region within the coronavirus RNA-dependent 
RNA polymerase (RdRp) gene (407 bp) accessed from GenBank, including a bat betacoronavirus sequence as an outgroup. Scale bar indicates the 
number of substitutions per site. Bootstrap support values are shown for nodes with 280% bootstrap support. Sequences from Vietnamese bats 


are indicated in red 
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KF294377|MI_bat_coronavirus|China 
KF294378/MI_bat_coronavirus|China 


400) KU836638!Porcine_epidemic_diarrhea_virus|United_Kingdom 
NC_003436!Porcine_epidemic_diarrhea_virus 
KX534206|Porcine_epidemic_diarrhea_virus|China 
KM077139!|Porcine_epidemic_diarrhea_virus|IUSA 
KU2526491Porcine_epidemic_diarrhea_virus|China 
KM609204|Porcine_epidemic_diarrhea_virus|China 
LC022792|Porcine_epidemic_diarrhea_virus|Japan 

100, KF294381IMd_bat_coronavirus!|China 
NC_028811|BatMr-AlphaCoVIChina 
97p 75_55_L13_R2_CoV_BAT_VN 
75_98 L13_R4 CoV_BAT VN 
75_55_L09_R3_CoV_BAT VN 
75 _55_L01_R4 CoV_BAT VN 
75 98 LO5 R6_CoV_BAT VN 
75_99 L01_R5_CoV_BAT_VN 
75_55_L03_R3_CoV_BAT_VN 
75_55 L10_R3_CoV_BAT_VN 
92] 75_55_L12_R6_CoV_BAT_VN 
AB539080/Bat_coronavirus|Philippines 
75_55_LO9_R6_CoV_BAT VN 
75_55_L13_R6_CoV_BAT VN 
75_55_L10_R6_CoV_BAT_VN 
75_98 L01_R3_CoV_BAT_ VN 
75_55_L15_R2 CoV BAT VN 
75_98 LO7_R2_CoV_ BAT VN 
75_98 LO5 R2 CoV BAT VN 
1009} 75 99 | CoV 
75_55_L01_R1_CoV_BAT_VN 
75_55_L06_R2 CoV_BAT_VN 
75_98 LO7_R4 CoV_BAT_VN 
75 55 _L03_R6_CoV_BAT VN 
75_55_L05 R5 CoV_BAT_ VN 
75_55 L15 R3_CoV_BAT VN 
75_55 LO04_R1_CoV_BAT_VN 
75_99 L12_R3_CoV_BAT_VN 
75_98 L15_R4 CoV_BAT_VN 
go] | 75_55_L03_R5_CoV_BAT VN 
75_99 LO2_R3_CoV_BAT_VN 
75_55_L02_R1_CoV_BAT_VN 
75_55_L04_R5 CoV_BAT_VN 
NC_009657|Scotophilus_bat_coronavirus_512|China 
DQ648821|Bat_coronavirus|China 
DQ648824|Bat_coronavirus|China 
75_55_L15_R5 CoV_BAT_VN 


100 


Lo6_R4_CoV_BAT_VN 


75_55_L11_R6_CoV_BAT VN 


75_55 L02_ R2 CoV_BAT VN 
75_55_LO5_R6 CoV_BAT_VN 
75_99 LO9_R1_CoV_BAT VN 
75_99 L15 R5 CoV_BAT VN 
75_99 LO3_R1_CoV_BAT_VN 
75_98_L14_R4 CoV_BAT_VN 
75 55 LO9 R5 CoV_BAT VN 
75_99 LO6_R1_CoV_BAT_VN 
NC_022103/Bat_coronavirusIUSA 
KF294382IMs_bat_coronavirus|China 
96 00 KJ4738 10)BatMs-AlphaCoV China 
99 NC_028814|BatRf-AlphaCoVIChina 
99 NC_018871|Rousettus_bat_coronavirus|China 
JQ989271|Rousettus_bat_coronavirus|China 
JQ989272!Hipposideros_bat_coronavirus_HKU10IHong_Kong 
JQ989266|Hipposideros_bat_coronavirus_HKU10|Hong_Kong 
JQ989269|Hipposideros_bat_coronavirus_HKU10|Hong_Kong 
HQ728481|Chaerephon_bat_coronavirus|Kenya 
A NC_028824|BatRf-AlphaCoViChina 
FE 203065|Rhinolophus_bat_coronavirus_HKU2 


_009988|Rhinolophus_bat_coronavirus_HKU2 
NC_028833/BatNv-AlphaCoVIChina 


100) DQ445911|Human_coronavirus_NL63|Netherlands 
KU521535|Human_coronavirus_NL63|IUSA 
AY518894|Human_group_1_coronavirus_associated_with_pneumonialNetherlands 
DJ009246|Human_coronavirus_NL63 
JX5241711IHuman_coronavirus_NL63|China 
KT253271|229E-related_bat_coronavirusiGhana 
100 KT2532701229E-related_bat_coronavirus|Ghana 
KT2532691229E-related_bat_coronavirus|Ghana 
KT253272|229E-related_bat_coronavirus|Ghana 
HOOg AF304460IHuman_coronavirus_229E 
X69721|Human_coronavirus_229E 
JX503060|Human_coronavirus_229E|Netherlands 
JX503061|Human_coronavirus_229Elltaly 
9 KT368908|Camel_alphacoronavirus|Saudi_Arabia 
JQ410000!Alpaca_respiratory_coronavirus|USA 
JIB KU291449|Camel_alphacoronavirus_229E|Kenya 
KT368894ICamel_alphacoronavirus|Saudi_Arabia 
KT368915I|Camel_alphacoronavirus|Saudi_Arabia 
41009 KJ473797|BatMf-AlphaCoVIChina 


NC_010438|Miniopterus_bat_coronavirus_HKU8IHong_Kong 
KJ473798|BatMf-AlphaCoVIChina 
100 KJ473799|BatMf-AlphaCoVChina 
KJ473800|BatMf-AlphaCoVIChina 
100 HQ728485IMiniopterus_bat_coronavirus|Kenya 
00 HQ728484IMiniopterus_bat_coronavirus|Kenya 
EU420137|Bat_coronavirus_1BlHong_Kong 
KJ473795|BatMf-AlphaCoV|China 
KJ473796|BatMf-AlphaCoV|China 
KJ473801 Baty pheCov i hing 
KJ473803/BatMf-AlphaCoV|China 
NC_010437/Bat_coronavirus_1AlHong_Kong 
KF294380!Rat_coronavirus|China 
NC_023760/Mink_coronavirus_strainIUSA 
0, LC119077|Ferret_coronavirus|Japan 
100 NC_030292/Ferret_coronavirus|Netherlands 
98 FJ938058|Feline_coronavirus|Netherlands 
GU553361|Feline_coronavirus|Netherlands 
99 HQ012369|Feline_coronavirus|Netherlands 
FJ938060!Feline_coronavirusIUSA 
KU215424|Feline_coronavirus|Belgium 
FJ938052|Feline_coronavirus|Netherlands 
FJ938057|Feline_coronavirus|Netherlands 
0.3 


EU186072IFeline_coronavirusIUSA 
substitutions/site 


FJ938053|Feline_coronavirus|Netherlands 
FJ938054|Feline_coronavirus|Netherlands 
FJ938056|Feline_coronavirus|Netherlands 
FJ938062|Feline_coronavirus|Netherlands 

KP143507|Feline_coronavirus|United_Kingdom 
DQ848678|Feline_coronavirus 
KP849472|Alphacoronavirus_1 Maly. 
GQ152141|Feline_coronavirus|Taiwan 
GQ477367\Canine_coronavirus|Taiwan 
700 76U856362!Canine_coronaviruslltaly 
KP981644|Canine_coronaviruslltaly 
04409|Canine_coronavirus 
75340ICanine_coronavirusIUSA 
JN634064/Feline_coronavirus 
JQ404410!Canine_coronavirus 
JQ408980!Feline_coronavirus 
KC175339|Canine_coronavirus|Germany 
JN856008|Canine_coronavirusIUSA 
KP202848!Transmissible_gastroenteritis_virus|China 
KU729220ITransmissible_gastroenteritis_virus|China 
FJ755618ITransmissible_gastroenteritis_virus|China 
DQ811787IPRCV_ISU-1IUSA 
KR270796!Porcine_respiratory_coronavirusIUSA 
LT545990|Swine_enteric_coronavirus 
NC_028806|Swine_enteric_coronaviruslltaly 
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HQ166910IZaria_bat_coronavirus|Nigeria 
NC_025217|Bat_Hp-betacoronavirus|China 


DQ084200!Bat_SARS_coronavirus 
400 KF294457|SARS-related_bat_coronavirus|China 


81 84 


94 99 


97 


84 
97 


100 


99 


0.2 
substitutions/site 


JX993987/Bat_coronavirus|China 
DQ648856IBat_coronavirus 
KJ47381 1|BatRf-BetaCoViChina 
KJ473812|BatRf-BetaCoViChina 
KJ473813|BatRf-BetaCoVIChina 
DQ648857!Bat_coronavirus 
KJ473814|BatRs-BetaCoVIChina 
JX993988IBat_coronavirus|China 
AY572035|Civet_SARS_coronavirus|China 
86 KF367457|Bat_SARS-like_coronavirus|China 
KT444582|SARS-like_coronavirus|China 
KJ473816|BatRs-BetaCoVIChina 
KJ473815|BatRs-BetaCoViChina 
FJ588686ISARS_coronavirus|China 
KF569996!Rhinolophus_affinis_coronavirus|China 
KC869678IBat_Coronavirus|South_Africa 
KU740200IMERS_coronaviruslEgypt 
KF917527IMERS_coronavirus|Saudi_Arabia 
KP719933IMERS_coronavirusIUAE 
KJ473821|BatVs-BetaCoV_SC2013!China 
100 KC545383|Hedgehog_BetacoronavirusiGermany 
NC_022643|Hedgehog_Betacoronavirus|Germany 
4100 _ KJ473822IBatTp-BetaCoViChina 
a5 DQ648794|Bat_coronavirus 


NC_009019!Tylonycteris_bat_coronavirus|China 
4100 g KJ473820IBatPa-BetaCoVIChina 


EF065510|Bat_coronavirus|China 
NC_009020!Pipistrellus_bat_coronavirus|China 
HQ728482I|Eidolon_bat_coronavirus|Kenya 
99 HM211098|Bat_coronavirus|China 
99 EF065516|Bat_coronavirus|China 
HQ728483|Rousettus_bat_coronavirus|Kenya 
DQ415911|Human_coronavirusIChina 
AY597011/Human_coronavirus 
KF686344|Human_coronavirusIUSA 
DQ415897!Human_coronavirus|China 
DQ415899|Human_coronavirus|China 
75_62_L12 R2_RAT_CoV_VN 
81 fF 75_63_Lo4_R5_RAT_CoV_VN 
75_65_L03_R6_RAT_CoV_VN 
75_65_L09_R2_RAT_CoV_VN 
75_63_L11_R6_RAT_CoV_VN 
400 KF294372|Longquan_RI_rat_coronavirus|China 
KF294370|Longquan_RI_rat_coronavirus|China 
KF294371 |Longquan_RI_rat_coronavirus|China 
AC_000192/Murine_hepatitis_virus_strain 
FJ647222|Murine_coronavirus|USA 
AF208067|Murine_hepatitis_virus 
FJ647220|Murine_coronavirus|USA 
91 FJ938068!Rat_coronavirus 
KF850449/Rat_coronavirusIUSA 
JQ173883|Murine_hepatitis_virus 
AB551247|Murine_hepatitis_virus|Australia 
FJ647223|Murine_coronavirusIUSA 
FJ647224|Murine_coronavirus|IUSA 
75_65_L02_R6_RAT_CoV_VN 
400 [P 75-85-L07_Ré_RAT_CoV_VN 


100 


99 


94 


KF294357ILongquan_Aa_mouse_coronavirus|China 


NC_026011|Rat_BetacoronavirusIChina 
75_62_L01_R6_RAT_CoV_VN 
75_65_L10_R5_RAT_CoV_VN 

100 _ JN874562/Rabbit_coronavirus|China 
JN874561|Rabbit_coronavirusIChina 
100 NC_017083|Rabbit_coronavirus|China 
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OUTGROUPIKJ473795|BatMf-AlphaCoVIChina 


DQ011855|Porcine_hemagglutinating_encephalomyelitis_virus|Belgium 


KF530093|Human_coronavirusIUSA 
KP198610IHuman_coronavirus|China 


JX860640ICanine_respiratory_coronavirus|South_Korea 


U00735|Bovine_coronavirus 
AF220295|Bovine_coronavirus 
LC061272|Equine_coronavirus|Japan 
EF446615|Equine_coronavirus 
LC061274|Equine_coronavirus|Japan 
AF391542IBovine_coronavirus 
FJ425190|Sambar_deer_coronavirus|USA 
FJ425187|White-tailed_deer_coronavirusIUSA 
FJ938065|Bovine_respiratory_coronavirus|USA 
FJ938063|Bovine_coronavirusIUSA 
KU558922|Buffalo_Betacoronavirus|Bangladesh 
FJ415324lHuman_enteric_coronavirus|Germany 
KF906249|Dromedary_camel_coronavirusIUAE 
EF424619|Bovine_coronavirusIUSA 
EF424621|Sable_antelope_coronavirus|IUSA 
EF424615|Bovine_coronavirusIUSA 
EF424623|Giraffe_coronavirusIUSA 


FIGURE 3. The phylogenetic relationships of betacoronaviruses sampled from Vietnamese rats. Maximum-likelihood phylogeny constructed 
using RNA sequences from nine Vietnamese rats and 81 reference sequences from the conserved region within the coronavirus RNA-dependent 
RNA polymerase (RdRp) gene (407 bp) accessed from GenBank, including a bat alphacoronavirus as an outgroup. Scale bar indicates the number 
of substitutions per site. Bootstrap support values are shown for nodes with 280% bootstrap support. Sequences from Vietnamese rats are 


indicated in blue 


(2016) recently observed similar characteristics for coronaviruses 
in Chinese bats. Phylogenetic analysis revealed that the sequences 


generated here belong to a recently detected and unclassified bat 


betacoronavirus lineage, and are closely related to other bat coro- 


naviruses sampled in China and the Philippines (Tang et al., 2006; 
Watanabe et al., 2010). 
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Using our sampling framework, we found that our ability to detect 
paramyxoviruses and coronaviruses was not associated with location 
or time of sample collection, and appeared to be both random and 
sporadic. Our findings indicate that the novel paramyxoviruses de- 
tected were not present in all bat colonies at the time of sampling, 
suggesting that the extent of carriage and/or excretion in any particu- 
lar bat colony is likely to fluctuate over time. These data indicated that 
while infection in a bat colony is periodic and transient for paramyxo- 
viruses, this may not be the case for coronaviruses. Therefore, from a 
zoonotic infection risk management perspective, it is prudent to as- 
sume that any bat colony could be infected at any time and that those 
working in proximity to the bat colonies should adopt exposure risk 
minimization strategies in case these viruses are capable of infect- 
ing humans. We additionally found that PCR amplification-positive 
samples from bats were associated with sample type; urine was 
more likely to be positive for paramyxoviruses while faecal material 
was more likely to generate positive amplification of coronaviruses. 
Importantly, the screening of these samples by PCR alone may under- 
estimate the prevalence of these viruses in the sampled populations, 
as detection is limited by the specificity of primers and the quality of 
samples collected using indirect methods such as that used for the 
collection of bat samples. 

Although none of the bat samples tested in this study indicate the 
presence of filoviruses, this does not provide conclusive evidence that 
these viruses are not circulating in Viet Nam. While fruit bats have been 
implicated as the primary reservoirs of filoviruses, a previous survey of 
multiple bat species in the Philippines confirmed Ebola Reston virus 
infection in one insectivorous bat species (Miniopterus schreibersii) of 
21 insectivorous and fruit-eating bat species sampled (Jayme et al., 
2015; Leendertz, Gogarten, Diix, Calvignac-Spencer, & Leendertz, 
2015). The samples included in the current study were collected only 
from the insectivorous S. kuhlii, while at least 90 other chiropterans 
have also been identified in Viet Nam (http://www. iucnredlist.org/). 
Our sampling was further limited by the use of faecal and urine sam- 
ples, which are unlikely to be optimum materials in which to detect 
filoviruses. Thus, a much larger and more comprehensive sampling of 
multiple bat species using more direct sampling methods would be 
necessary to detect and understand the potential circulation of filo- 
viruses in the region. 

In contrast to previous investigations (Woo et al., 2012), paramyxo- 
virus nucleic acid was not detected in any the rat faecal samples 
tested. A possible explanation as to why the rats tested were all nega- 
tive for paramyxoviruses could be that the sample type (faeces) is not 
appropriate to detect the virus. This observation is supported by a pre- 
vious study, where a novel rat paramyxovirus was detected only in the 
kidney and spleen of rats (Woo et al., 2012). Novel coronavirus RNA 
was detected by PCR amplification and characterized by sequencing 
in 4.4% of the rat faeces. The obtained rat coronavirus sequences 
were closely related to rat and murine betacoronavirus sequences 
from southern China that were found to be capable of infecting up 
to ten different rat species, including R. losea and R. tanezumi (Wang 
et al., 2015). In our study, novel rat coronavirus strains were detected 


in a single rat species (R. argentiventer). These data do not exclude the 
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possibility that other rat species carry this coronavirus virus, particu- 
larly as our sampling was biased towards R. argentiventer (234/270). 

Despite a lack of appropriate control measures for zoonotic vi- 
ruses in Viet Nam, viral spillover from bats and/or rats has not, as yet, 
been observed, although these preliminary findings suggest that both 
coronaviruses and paramyxoviruses may be endemic in Vietnamese 
bat populations. While we are unable to comment on the pathogenic 
potential of these novel viruses in humans, it is likely that if such infec- 
tions do arise, they are rare and go undiagnosed. Paramyxoviruses in 
particular are known to cause infections in humans, and the spillover 
of Hendra virus and Nipah virus from bats to horses or pigs, respec- 
tively, and subsequently to humans has been reported in Australia 
and Malaysia (Edson et al., 2015; Goldspink et al., 2015; Sherrini & 
Chong, 2014). Laboratory investigations are being conducted to as- 
sess whether the coronaviruses and paramyxoviruses detected here 
are viable in cell culture in order to better predict whether these vi- 
ruses have the potential to cause spillover infections in the human 
population and to improve genomic characterization of these viruses 
to allow for greater phylogenetic inference. Furthermore, our country- 
wide project is aiming to recruit hospitalized patients with common 
and uncommon disease phenotypes in order to determine whether 
such viruses cause previously undiagnosed infections in humans in 
Viet Nam (Rabaa et al., 2015). 

In conclusion, we have detected the presence of bat paramyxovi- 
ruses, bat coronaviruses and rodent-associated coronaviruses within 
a limited geographic area and time frame in Viet Nam. Consequently, 
we suggest that larger-scale surveillance is needed to fully understand 
the role played by rodents and bats in the evolution, emergence and 


dissemination of paramyxoviruses and coronaviruses in Viet Nam. 
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